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****************************************************************************************
DISCLAIMER
This book contains a number of computer programs written in Matlab version 5.2.  
They are intended to use by the readers to understand the concepts of adsorption better, 
and they are for educational and instructional purposes only.  Use of these programs for 
design purposes may be done at the owner's risk.  Neither the author nor the publisher 
will bear any responsibility for the accuracy obtained from the enclosed programs.
****************************************************************************************

Program 	Program description
-------		-------------------
ADSORB0		Adsorption kinetics of one species in a particle
			1. Linear isotherm
			2. Isothermal conditions
ADSORB1A	Adsorption kinetics of one species in a particle
			1. Langmuir
			2. Isothermal conditions
ADSORB1B	Adsorption kinetics of one species in a particle
			1. Langmuir isotherm
			2. Nonisothermal conditions
ADSORB1C	Adsorption kinetics of one species in a crystal
			1. Langmuir isotherm
			2. Isothermal conditions
ADSORB1D	Adsorption kinetics of one species in a crystal
			1. Langmuir isotherm
			2. Nonisothermal conditions
ADSORB1E	Adsorption kinetics of one species in a zeolite pellet composed of many crystals
			1. Langmuir isotherm
			2. Nonisothermal conditions
			3. Zeolite pellet = many  crystals
ADSORB3A	Adsorption kinetics of N species in a crystal under isothermal conditions
			1. Extended Langmuir equation
			2. Isothermal conditions
			3. N species
			4. Crystal
ADSORB5A	Adsorption kinetics of N species in a cylindrical pore under isothermal conditions
			1. Extended Langmuir equation
			2. Isothermal conditions
			3. Cylindrical pore
			4. N species
			5. Viscous flow
ADSORB5B	Adsorption kinetics of N species in a particle under isothermal conditions
			1. Extended Langmuir
			2. Isothermal conditions
			3. Porous particle
			4. N species
			5. Viscous flow and surface flow
ADSORB5C	Adsorption kinetics of N species in a particle under non-isothermal conditions
			1. Extended Langmuir
			2. Non-isothermal conditions
			3. Particle
			4. N species
			5. Including viscous flow and surface flow
BIND		Calculation of binary Dij
CAPILL2		Steady state diffusion flux in a capillary
			1. Knudsen flow
			2. Molecular diffusion flow
E_VS_D		Interaction energy versus half-width or pore radius
EIGENP		Eigenvalue in a single particle
FASTIAS		Fast IAS theory
FOWLER		Adsorption isotherm of Fowler-Guggenhein equation
HILL		Adsorption isotherm of Hill-de Boer equation
HK		Effective width of a slit pore as a function of reduced pressure using HK method
IAS		Ideal Adsorption Theory
ISO_FIT1	Pure component isotherm fitting
ISO_FIT2	Pure component fitting at multiple temperatures
MesoPSD1	Mesopore size distribution using double gamma PSD
P1043		Characteristics of 10_4_3 potential energy
PSD_MESO	Statistical thickness as a function of reduced pressure for mesopore
STEFTUBE	Stefan tube of a mixture of n species evaporating into an n-th species medium
UPTAKEP		Uptake calculation in a single particle


Most of these program can be run by typing the name of the program.  For programs 
which require optimization routines, the user will be prompted for the guessing parameters 
and for approval whether the initial guess is acceptable.  To change the parameters in the 
program, the user simply open the m files using either the NOTEPAD in Windows or the built-in 
m-editor of the MatLab version 5.2.  The built-in m-editor is recommended because it can help 
the user whether the syntax used is correct or not.  In each program there is a section called 
"USER SUPPLY SECTION".  There is a nomenclature in each of the program so that the user knows 
the meaning of each parameter.  The user may change these parameters and save the file before 
executing the program within the MatLab worksheet.

For example, if the user wants to simulate the adsorption kinetics of an adsorbate in a 
particle under isothermal conditions and Langmuir isotherm, the user simply opens the 
ADSORB1A.M file, and will see the following "USER SUPPLY SECTION"

 %---------------------------------------------------------------------------------------------------
% USER SUPPLY SECTION
s		= 0;			% particle shape factor
R		= 0.2;			% particle radius (cm)
porosity	= 0.33;			% particle porosity (-)
ci		= 0;			% initial gas concentration (mole/cc)
cb		= 1e-6;			% bulk gas concentration (mole/cc)
c0		= 1e-6;			% reference concentration (mole/cc)
b		= 1e6;			% adsorption affinity (cc/mole)
cmus		= 5e-3;			% maximum adsorption capacity (mole/cc)
biot		= 1000;			% biot number (-)
Dp		= 0.1;			% pore diffusivity (cm^2/sec)
Ds		= 1e-5;			% surface diffusivity (cm^2/sec)
%-----------------------------------------------------------------------------------------------------

To change any of those parameters, for example the particle radius to 0.1, the user simply 
replace R = 0.2 in the "USER SUPPLY SECTION" to R = 0.1.  Next save the file, and then return 
to the MatLab worksheet and execute the program by simply typing ADSORB1A.


SUPPORTING PROGRAMS
All the programs enclosed with this book require a number of supporting programs to run.  
Most of these programs are built-in programs within the MatLab version 5.2.  However, 
some programs such as the Newton-Raphson program for solving nonlinear algebraic equations 
and collocation programs for solving boundary value problems are enclosed with the CD.  
The collocation programs were adapted from algorithms developed by Villadsen and 
Michelsen (1978) (Chapter 3).  Readers must refer to this excellent treatise for further 
exposition of the orthogonal collocation method.  


BASIC INSTRUCTIONS
Prior to running the programs enclosed in the diskette, the user should copy them into a 
sub-directory in a hard drive and make sure that the paths of these programs are added to 
the MatLab worksheet so that they can be executed.
